SAIC-98/1010 


“Global Magnetohydrodynamic 
Modeling of the Solar Corona” 

FINAL REPORT 
( 2 / 07 / 95 — 2 / 06 / 98 ) 



Science Applications International Corporation 

An Employee-Owned Company 


NASA CONTRACT: NASW-4968 
Space Physics Supporting Research and Technology (SR&T), 
and Suborbital Program 


Principal Investigator: 

Jon A. Linker 

Science Applications International Corporation 
10260 Campus Point Drive 
San Diego, CA 92121-1578 


10260 Campus Point Drive , San Diego , California 92121 (619) 546-6000 

Other SAIC Offices: Albuquerque, Boston, Colorado Springs, Dayton, Huntsville, Las Vegas, Los Angeles, McLean, Oak Ridge, Orlando, Palo Alto, Seattle, and Tucson 


Global Magnetohydrodynamic Modeling of the Solar Corona 

NASA Contract NASW-4968 (2/07/95-2/06/98) 

P.I.: Jon Linker 
Final Report 

1. Introduction 

The coronal magnetic field defines the structure of the solar corona, the position 
of the heliospheric current sheet, the regions of fast and slow solar wind, and the 
most likely sites of coronal mass ejections. There are few measurements of the 
magnetic fields in the corona, but the line-of-sight component of the global magnetic 
fields in the photosphere have been routinely measured for many years (for example, 
at Stanford’s Wilcox Solar Observatory, and at the National Solar Observatory 
at Kitt Peak). The SOI/MDI instrument is now providing high-resolution full- 
disk magnetograms several times a day. Understanding the large-scale structure of 
the solar corona and inner heliosphere requires accurately mapping the measured 
photospheric magnetic field into the corona and outward. 

Ideally, a model should not only extrapolate the magnetic field, but should 
self- consistently reconstruct both the plasma and magnetic fields in the corona and 
solar wind. Support from our NASA SR&T contract has allowed us to develop 
three-dimensional magnetohydrodynamic (MHD) computations of the solar corona 
that incorporate observed photospheric magnetic fields into the boundary condi- 
tions. These calculations not only describe the magnetic field in the corona and 
interplanetary space, but also predict the plasma properties as well. Our compu- 
tations thus far have been successful in reproducing many aspects of both coronal 
and interplanetary data, including the structure of the streamer belt, the location 
of coronal hole boundaries, and the position and shape of the heliospheric current 
sheet. 

The most widely used technique for extrapolating the photospheric magnetic 
field into the corona and heliosphere are potential field models, such as the potential 
field source-surface model (PFSS) ( Schatien et al ., 1969; Altschuler and Newkirk , 
1969; Hoeksema ., 1984, 1991) and the potential field current-sheet (PFCS) model 
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( Schatien , 1971). 


The potential field models have been shown to describe many aspects of coro- 
nal and interplanetary data (e.g., Wang and Sheeley , 1988; Hoeksema and Suess , 
1990: Wang and Sheeley , 1992), and the models are still being improved ( Schulz , 
1995; Z/mo and Hoeksema , 1995). Ease of computation makes these models es- 
pecially useful tools, but the simplifying assumptions underlying these methods 
limit their applicability. The coronal phenomena of interest are large-scale and 
long-wavelength, and are thus well described by the magnetohydrodynamic (MHD) 
equations. Magnetostatic solutions ( Bogdan and Low , 1986; Bagenal and Gibson , 
1991; Gibson and Bagenal , 1995; Gibson et a/., 1996) can describe the density 
distribution in the lower corona but cannot model the solar wind. Full MHD com- 
putations (including flow terms) contain much of the physics required for describing 
the solar corona and inner heliosphere. At the inception of our NASA-supported 
project, MHD computations of helmet streamers had been performed for many years 
(e.g., Endler , 1971; Pneuman and Kopp , 1971; Steinolfson et aZ., 1982; Washimi et 
al, 1987; Linker et al, 1990; Wang et al., 1993; Linker and Mikic , 1995; Stewart 
and Bravo , 1996). However, these calculations typically used idealized magnetic 
flux distributions (such as a dipole or other multipole component), and most were 
limited to two dimensions. The resulting configurations could only be compared to 
the corona in terms of abstract qualities, and not specific observations. 

To perform a realistic 3-D MHD computation of the corona, it is necessary 
to incorporate realistic magnetic fields into the boundary conditions ( Mikic and 
Linker , 1996; Linker et al. , 1996; Usmanov , 1996). We have developed numerical 
MHD computations of the corona and inner heliosphere that incorporate observed 
photospheric magnetic fields. In Section 2, we briefly describe the methodology of 
our computations. Section 3 shows comparisons of our results with specific obser- 
vations. Section 4 describes improvements to the model that are presently being 
investigated. 

2. Modeling the Corona and Inner Heliosphere 

To compute self-consistent three-dimensional MHD solutions for the large-scale 
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corona, we solve the following form of the equations in spherical coordinates (Mikic 
and Linker , 1994. describes the method): 
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where B is the magnetic field intensity, J is the electric current density, E is the 
electric field, v, p, and p are the plasma velocity, mass density, and pressure. The 
gravitational acceleration is g, 7 is the ratio of specific heats, 77 is the resistivity, v is 
the viscosity, S represents energy source terms, and the wave pressure p w represents 
the acceleration due to Alfven waves [ Jacques , 1977; Hollweg , 1978]. 

The term S in equation (6) includes the effects of coronal heating, thermal 
conduction parallel to B, radiative losses, and Alfven wave dissipation (viscous and 
resistive dissipation can also be included). 5 and P w were set to zero for the results 
discussed in this section, which restricts us to examining polytropic solutions (i.e., 
dPp-~ / dt — 0). These solutions have the advantage that relatively simple models 
can match many of the properties of the corona: however, values of 7 close to 1 
(7 — 1.05 for the results shown here) are necessary to produce plasma profiles that 
are similar to coronal observations ( Parker , 1963). This reflects the importance of 
the processes embodied in 5 that we have neglected here. Section 4 describes the 
incorporation of these processes. 

The methods we use to solve equations (1-6), including the boundary condi- 
tions. have been described previously ( Mikic and Linker , 1994, 1996; Linker et al. 
1996: Linker and Mikic , 1997). We have used synoptic magnetograms from the 
Wilcox Solar Observatory at Stanford and the National Solar Observatory at Kitt 
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Peak (collected during a solar rotation by daily measurements of the line-of-sight 
magnetic field at central meridian) to specify the radial magnetic field (B r ) at the 
photosphere (in the manner described by Wang and Sheeley , 1992). We also spec- 
ify the density and pressure at the lower boundary. For the initial condition, we 
compute the potential (current-free) field consistent with the specified distribution 
of B r at the lower boundary, and a wind solution [Parker, 1963] consistent with the 
specified p and p. We then solve equations (1-6). The configuration evolves until a 
steady state is reached, when the magnetic field lines and plasma state have settled 
into equilibrium. A typical run of our code on a relatively coarse (nonuniform) com- 
putational mesh (81 x 51 x 32 r,6U mesh points) takes about an hour of Cray-C90 
CPF time per day of real time, while higher-resolution runs (101 x 101 x 64) take 
a few hours of CPU time per day of real time. Coronal computations relax to a 
steady state in a few days of real time. 

3. Comparisons with Observations 

The solutions obtained in the manner described in Section 2 can in principle 
provide a 3-D description of the corona and inner heliosphere, including the detailed 
distribution of magnetic fields, currents, plasma density, and temperature. The 
validity of this approach can only be verified through comparison with observations. 
With this goal in mind we have sought to compute solutions for specific time periods 
of interest. Carrington rotation 1S92 (CR1892, January 27-February 23, 1995) is 
good example: a variety of coronal data is available for that time period (Mauna Loa 
white light images. \ ohkoh soft X-ray images, and coronal hole boundaries based 
on the Kitt Peak 10830 images). At the same time Ulysses probed the heliosphere 
over a wide range of latitudes. We used a Kitt Peak synoptic map to specify B r in 
the manner described in Section 2, and computed a solution out to 400 solar radii 
(1.9 astronomical units, or A.U.). Figure 1 shows how the solution captures coronal 
structures near the Sun and also models the inner heliospere. 

Coronal or helmet streamers appear as bright regions in white-light coronal im- 
ages; they are believed to outline closed magnetic regions on the Sun. To directly 
compare our results with observations, we must develop images like those obtained 
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Magnetohydrodynamic Modeling of the Solar Corona and Inner Heliosphere 
Carrington Rotation 1892 (January 27 - February 23, 1995) 
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Figure 1. An MHD calculation of the solar corona and inner heliosphere for Carrington Rotation CR1892, showing (a) magnetic field lines (IBI is 
contoured on the solar surface) (b) the polarization brightness (pB) (c) mesh close to the Sun (d) boundary between open field regions (coronal holes, 
shown as dark gray) and closed field regions (light gray) (e) spiral field lines viewed from 1.9 A.U. above the pole of the Sun (f) mesh at 1.9 A.U., with an 
isosurface showing the position of the heliospheric current sheet. The variety of coronal and interplanetary properties predicted by the calculation allow 
the results to be tested against a wide range of data sets. 





with coronagraphs and during eclipses. Frequently the polarization brightness (pB) 
is the observed quantity. The distribution of pB in the plane of the sky is pro- 
portional to the line-of-sight integral of the product of the electron density and a 
scattering function that varies along the line of sight ( Billings , 1966). Using the 
plasma density from our computation, we can calculate pB and compare it with 
actual observations. Eclipse and coronagraph images typically compensate for the 
rapid fall off of coronal density, through the use of radially graded filters or the 
vignetting properties of the coronagraph. We detrend our computed pB in a similar 
manner. Figure 2 shows a comparison between the predicted pB from our computa- 
tion with Mauna Loa Coronameter data. Throughout the rotation, the results from 
the MHD computation correspond well with the observations. This computation, 
along with those we have performed for other time periods, support the long-held 
belief that the magnetic field distribution on the Sun controls the position and shape 
of the streamer belt. 

Open field regions on the Sun (coronal holes) are the source of the solar wind. 
Another way we can test the MHD model is to compare open field regions from 
the computation with coronal holes observed from solar images. By tracing field 
lines at each grid point and determining whether the field lines were open (reached 
the outer boundary) or closed (returned to the Sun), we mapped the open field 
regions predicted by our computation and compared them with Yohkoh soft X-ray 
images, as shown in Figure 3. Coronal holes typically show up as dark regions 
in X-rays. Throughout the rotation, the coronal hole boundaries predicted by the 
model agree reasonably well with the dark regions seen in Yohkoh, although there 
are also some discrepancies. It is important to note that the identification of coronal 
holes in soft X-ray images (or the Kitt Peak 10830 line) is somewhat subjective. 
In X-rays, neighboring coronal emission can obscure coronal holes, and the absence 
of emission can also occur in quiet closed field regions. As part of a coordinated 
study for the time period of the Ulysses fast latitude scan {Neugebauer et al . , 1997), 
the coronal hole boundaries from our model were also compared with the He 10830 
observations from Kitt Peak. Reasonable correspondence was again found, although 
it is interesting to note that not only were there some disagreements between the 
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Comparison of a 3D MHD Model of CR1892 with Mauna Loa Coronameter Data (HAO) 

January 28, 1995 February 4, 1995 February 10, 1995 February 14, 1995 February 21, 1995 



Observations of Polarization Brightness (Mauna Loa MK3) 



Polarization Brightness (MHD Model) 



Field Lines (MHD Model) 


Figure 2. Comparison of Mauna Loa Coronameter data (top panels) with the polarization brightness (pb) computed from die MHD 
simulation (middle panels). The model results correspond well to the observations throughout the rotation. In the bottom set of panels 
traces of the magnetic field lines are shown, with the magnitude of the magnetic field contoured on the solar surface. The strongest 
fields (yellow colors) correspond to active regions. 


Comparison of Coronal Hole Boundaries: A 3D MHD Model of CR1892 and Yohkoh Data 


January 28, 1995 February 4, 1995 February 11, 1995 February 15, 1995 February 22, 1995 



Yohkoh Soft X-Ray Data 



Field Lines and Coronal Hole Boundaries (MHD Model) 


Figure 3. Comparison of Yohkoh soft X-ray images (top panels) with the open field regions predicted by the MHD simulation. 
Identification of coronal hole boundaries in Yohkoh data can be somewhat subjective. The darkest regions often correspond to coronal 
holes (although the absence of emission can occur in closed field regions, and coronal holes can be obscured by neighboring emission). 
Nevertheless, the coronal hole boundaries predicted by the MHD calculation (shown in dark gray) correspond approximately to those 
seen in the Yohkoh data. 
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coronal holes predicted by the model and the absence of emission in the observations, 
coronal hole boundaries in He 10830 and the dark regions in Yohkoh also did not 
always correspond. 

The comparisons of our results with coronal data, shows that the basic large-scale 
structure of the corona has been captured in the model. We have also compared 
our results with interplanetary data. During February-April, 1995 (Carrington ro- 
tations 1892-1894), the Ulysses spacecraft sampled a wide range of heliographic 
latitude in a relatively short period of time; this time period has since been referred 
to as the fast latitude scan. Figure 4 (red curve) shows “simulated” Ulysses data, 
created by flying the Ulysses trajectory through the MHD solution and extracting 
D, along the trajectory. For comparison, the actual B T measured by Ulysses is also 
shown (blue curve). The fluctuations in the solar wind data are caused by turbulent 
waves and transient phenomena, and thus are not present in the MHD model, but 
the model does reproduce the correct magnitude and polarity of the field, includ- 
ing many of the heliospheric current sheet crossings. Not surprisingly, the model 
matches the data best during and shortly after CR1892, when the photospheric 
data used to perform the calculation is most applicable. Photospheric data for the 
subsequent rotations shows that the magnetic field is indeed evolving. 

Figure 5 shows the heliospheric current sheet predicted by our computation for 
CR1892. depicted as an isosurface. The Ulysses trajectory in the rotating frame of 
the Sun is shown. Blue colors for the trajectory correspond to times when negative 
B r fields were observed, and red colors show times when positive B r fields were 
observed. The dominant polarity observed by Ulysses agrees with that predicted 
bv the simulation for most of the early (lower) part of the trajectory. In the later 
part of the trajectory, there is a passage into negative polarity not predicted by the 
model (this can also be seen in Figure 4 around day-of-year 85). As this occurred in 
C'R 1894. it is possible that the heliospheric magnetic field has changed significantly 
from C'R 1892 (two months earlier). To test this possibility, we computed MHD 
models for rotations CR1893 and CR1894 as well. Figure 6 shows a comparison 
of heliospheric current sheet (HC’S) crossings observed by Ulysses and projected 
back to the Sun (Smith et al. 1995) with the HCS predicted by the MHD model for 


8 



Magnetic Field Comparison 
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Figure 4. A comparison of a simulated time series of B r for the MHD computation of CR1892 (red curve) with Ulysses data 

(blue curve). Heliospheric current sheet crossings identified by Smith et al (1995) are marked as crosses (transient phenomena 
can also cause B r to change sign). The time period for CR1892 (when the photospheric data used by the MHD calculation is 

most applicable) is marked on the plot. The MHD calculation reproduces the correct magnitude of the field and matches many 
of the current sheet crossings. Not surprisingly, the MHD results match the data best during or close to the CR1892 time 
period. 



Ulysses Trajectory Shown with a Heliospheric 
Current Sheet Computed Using a 3D MHD Model 



Heliospheric Current Sheet 

Fast Latitude Scan (Jan. - Apr. 1995) 

Carrington Rotations 1892, 1893, 1894 

Figure 5. The heliospheric current sheet (HCS) predicted by the MHD computation of CR1892 shown in Figure 1. The 
Ulysses trajectory (red and blue curve) is shown in the rotating frame of the Sun. The colors of the trajectory indicate the 
polarity of the magnetic field observed by Ulysses. Fluctuations in B r are present because of waves and transients in the solar 

wind, but the dominant polarity observed by Ulysses is similar to that predicted by the MHD computation. 


Comparison of MHD Computations of CR1892, CR1893, and CR1894 
with Ulysses Heliospheric Current Sheet (HCS) Crossings 



Figure 6. The HCS predicted by MHD computations of CR1892, CR1893, and CR1894 (red, blue, and green curves) and the 
Ulysses trajectory during this time (red, blue, and green lines) projected back at the Sun. HCS crossings from Ulysses 
identified by Smith et al (1995) (also projected back to the Sun) are shown as crosses. There is overall good agreement 
between the predicted HCS crossings (where the Ulysses trajectory intersects the computed HCS) and the observed crossings. 
The MHD calculations predict that the HCS moved higher in latitude in the north during CR1894, as seen in the data. 
However, the uppermost crossings observed by Ulysses are not predicted by the model. Sensitivity of the calculations to the 
polar fields (poorly determined by line-of-sight magnetograms) could be a possible reason for this difference. 
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the two rotations. The red, blue and green straight lines show the Ulysses trajec- 
tory for CR1892, CR1893, and CR1894, respectively, and the red, blue, and green 
curves show the predicted HCS for these rotations. There is overall good agreement 
between the predicted HCS crossings (where the Ulysses trajectory intersects the 
model HCS) and the observed crossings (marked as +). The MHD calculations 
based on the photospheric data for CR1893 and CR1894 predict that the helio- 
spheric magnetic field should have evolved during this time period and that the 
HCS should have moved to higher northern latitudes for CR1894, consistent with 
the observations. However, the two uppermost crossings are still not predicted by 
the model. We note that the tilt of the HCS (and thus its extent in northern and 
southern latitudes) can be strongly influenced by the polar fields, which are not well 
determined by line-of-sight magnetograms. 

We have also computed MHD solutions for a number of other time periods 
of interest, including the last three total solar eclipses (see Figure 7). These and 
other comparisons of our MHD solutions with observations show that the MHD 
model typically does a good job of reproducing many aspects of coronal structure. 
While the success of the model has been encouraging, there is still significant room 
for improvement. In the next section, we describe how we are working refine and 
extend the physics in our computations to provide a tool capable of addressing 
important aspects of coronal structure that are beyond the scope of the present 
model. 

4. An Improved Model of the Corona and Solar Wind 

Starting at the photosphere and rising upward into the solar atmosphere, the 
temperature rises steeply in the chromosphere and transition region. The detailed 
coronal heating mechanism is not yet understood (see Parker , 1994 for a recent 
review), but most coronal heating mechanisms imply that this sharp temperature 
rise is the result of a volumetric heat source in the corona. In the inner corona, the 
large parallel thermal conductivity tends to make the temperature relatively uniform 
(on the order of 1-2 xlO 6 °K). The density at the top of the transition region is 
determined by the balance between radiation loss in the chromosphere and heating 
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March 9, 1997 


LASCO C2 Coronagraph Image 


Figure 7. Comparison of MHD computations of the solar corona with eclipse observations. For the November 3, 1994 solar eclipse, we used a 
synoptic magnetogram for Carrington rotation 1888 (October 10 - November 6, 1994) as boundaiy conditions for the model, and we compared 
with the results subsequent to the eclipse. For the October 24, 1995 and March 9, 1997 solar eclipses, we used synoptic magnetograms from the 
rotation prior to the eclipse (CR1900: September 2 - September 29, 1995, and a combination of CR1918 and 1919: January 30 - February 26, 
1997) to predict the structure in advance. All of these calculations used synoptic magnetograms from the Wilcox Solar Obsevatory . For the March 
9 eclipse, no eclipse data is yet available. Comparison with an image from the LASCO coronagraph aboard the SOHO spacecraft is shown. These 
and other coronal modeling results can be found on our WEB page: http://iris023. saic.com: 8000/corona/modeling.html 





by thermal conduction from the hot corona ( Withbroe, 1988). As we extend into the 
outer solar corona and interplanetary medium, the temperature decreases slowly as 
a result of solar wind acceleration and thermal conduction losses. Beyond about 
10 I?., (where R s is the solar radius), the plasma becomes collisionless, reducing 
the thermal conduction (Hollweg, 1978). In this region, solar wind acceleration by 
Alfven waves can be important, and may be necessary to produce the observed fast 
solar wind at 1 A.U. 

For the solutions described in Section 2, the use of a polytropic energy equation 
with a reduced polytropic index 7 is an attempt to combine all these effects into 
a simplified energy equation. However, not surprisingly, this simple model fails to 
reproduce the fast (« 800 km/s) and slow (ss 400 km/s) wind streams that are 
measured in interplanetary space. The nearly uniform coronal temperature that 
results from this model implies that the density scale height does not change much 
between the open and closed field regions, and so the magnitude of the density 
contrast between coronal holes and streamers is also not reproduced. 

It is the complicated interplay of radiation loss, thermal conduction, coronal 
heating, and Alfven wave acceleration and dissipation that describes the accelera- 
tion of the solar wind from the inner solar corona into the heliosphere. Fortunately, 
we can benefit from the considerable amount of research that has already been per- 
formed with one-dimensional MHD and multi-fluid models ( Withbroe, 1988; Habbal 
et ai. 1995). These models have been quite successful, despite their obvious geo- 
metrical limitations, in modeling thermodynamic processes in the solar wind and 
in comparisons with spacecraft solar wind measurements (Esser and Habbal , 1995; 
Habbal et al . , 1996). The prescription for improving our coronal model follows nat- 
urally from existing 1 -D models, and requires the inclusion of the effects described 
above in the energy and momentum equations of the 3 -D MHD model. Some of 
these effects (heating and thermal conduction) have previously been included in 
2 -D MHD models (Suess et al.. 1996). The advantage of a multi-dimensional MHD 
ealailation is that the magnetic field geometry (i.e., the location and distribution 
of open and closed field regions) is determined self-consistently, eliminating the 
parameters in the 1 -D models that specify the flux tube geometry (the so-called 
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“expansion factor” ). As has been done in 1-D models, the inner boundary is now 
chosen to be within the transition region, at a temperature of To = 500, 00Q°K. 
The balance of radiation loss and thermal conduction within the chromosphere and 
transition region determines the density at the inner boundary from the condition 
( Withbroe , 1988) 

Cn-r 0 2 [ ° Q{T)T 1/2 dT = i, 0 2 (To) (7) 

JT c h 2 

where n e is the electron density at the inner boundary, C is a known constant, 
<70 (To) is the thermal conduction heat flux at To, and the integral is performed from 
the base of the chromosphere (at T c h % 6000° Iv) through the top of the transition 
region (at T = To). The boundary conditions on the velocity are determined from 
the characteristic equations ( Linker and Mikic , 1995, 1997: Mikic and Linker , 1996). 
Therefore, in this formulation, the only boundary conditions required from obser- 
vations at the lower boundary r = R s are on the radial magnetic field. 

The energy equation (6) becomes 

dp 

-T + V • (pv) = (7 - l)(-pV • v - V • q - n e n p Q(T) + H ch + D) (8) 

where q = — Kybb • VT is the parallel (along B) thermal conduction heat flux, 
Hch is the coronal heating source, D is the Alfven wave dissipation term, n e and 
n p are the electron and proton density, and Q(T) is the radiation loss function 
(Rosner et al ., 1978). The parallel thermal conductivity Ky is the Spitzer value 
in the collisional regime, fey = 9 x 10'T 5 / 2 [erg/cm 2 /s], with T in degrees K, and 
is reduced in the collisionless regime (beyond % 10 R s )\ the representation in this 
regime is parameterized as a. collisionless heat flux (Hollweg. 1978). 

Results from previous 1-D computations suggest that an additional source of 
energy and momentum is necessary to accelerate the solar wind. Alfven waves of 
solar origin are considered to be the most likely candidate ( Hollweg , 1978; Withbroe, 
1988). The acceleration of the solar wind by Alfven waves occurs on a spatial and 
time scale that is below the spatial and time resolution of a global numerical model. 
In 1-D models, this effect is included using an equation for the time-space averaged 

Alfven wave energy density e ( Jacques , 1977) 

de 

— + V • F = v • Vp u . - D (9) 
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where F = (|v + v A )e is the Alfven wave energy flux, v A = B/^/Wp is the Alfven 
speed, v A = is the Alfven velocity, and p w = is the Alfven wave pres- 

sure that appears in the momentum equation (5) and represents the force on the 
plasma by Alfven waves. To incorporate equation (9) in a multi-dimensional imple- 
mentation. it is necessary to transport two Alfven wave fields (waves parallel and 
antiparallel to B), which are combined to give the Alfven wave energy density e. e is 
related to the space-time average of the fluctuating component of the magnetic field 
SB by e = (SB~) /Air. The dissipation term D expresses the nonlinear dissipation 
of Alfven waves in interplanetary space and will be modeled phenomenologically 
( Hollweg , 1978). The term v • Vp w in eq. (9) is the work done on the Alfven waves 
by the plasma flow. 

The increased complexity of the improved model makes a step-by-step approach 
essential. We first incorporated our improved model in a 1-D code, and reproduced 
the results of Withbroe (1988). We have since implemented the model in two di- 
mensions, and performed a preliminary comparison between a polytropic model, 
a thermodynamic model (heating, parallel thermal conduction, and presence of a 
transition region), and a. thermodynamic model with Alfven waves. Figure 8 shows 
this comparison. With the improved thermodynamic model, the observed white 
light contrast between coronal holes and streamers is reproduced ( Bagenal et dl. 
1996), and the wind achieves faster velocities. With Alfven waves, a larger contrast 
between the fast and slow wind develops. 

The improved model significantly advances our capabilities to study the large- 
scale corona and inner heliosphere. The trade-off is that some of the important 
physical processes are poorly understood (coronal heating, Alfven wave dissipation, 
and collisionless thermal conduction) and must be parameterized. In the results 
shown in Figure 8. the coronal heating was parameterized as 

H ch = #o(0)exp[- (r - R s ) / A(0)] (10) 

where Ho{9) expresses the latitudinal variation of the volumetric heating, and \{6) 
expresses the latitudinal variation of the heating function scale length. For the cases 
shown in Figure S. Ho. A. and the flux of Alfven waves were assumed to be uniform 
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(a) 



Polytropic Model 
(y = 1.05) 


Thermodynamic Model Thermodynamic Model 
(no Alfven waves) (with Alfven waves) 



Polarization Brightness 




Figure 8. Pr elimi nary comparison between a polytropic model and a thermodynamic model (with and without Alfven waves) for 
a 2-D simulation with an initial dipole magnetic field, (a) The observed contrast in polarization brightness between coronal holes 
and streamers is reproduced by the models with improved thermodynamics, (b) With heating and thermal conduction present, the 
closed field region becomes hotter than the open field, (c) The thermodynamic models predict faster solar wind velocities, 
especially when Alfven wave acceleration is present. With Alfven waves, a greater contrast between the wind speed at high 
latitudes and the equator develops. 
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in latitude. In the future, we plan to study the effects of the different parameters in 
2-D, and change the coronal heating to be expressed in terms of magnetic topology 
(i.e.. a proxy for open and closed field regions) rather than latitude 9. Eventually 
we plan to incorporate the improved thermodynamics into full 3-D solutions, just 
as we have previously done with the polytropic model. 
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